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1.0 Introduction 

Multiuser communications systems that employ code division multiple 
access (CDMA) exhibit a user capacity limit in the sense that there exists a 
maximum number of users that can simultaneously communicate over the 
channel for a specified level of performance per user. This limitation is 
brought about by the ultimate domination of the other user mterference over 
the additive thermal noise. Over the years researchers have sought ways to 
extend the user capacity of CDMA systems either by employing optimum 
[maximum-likelihood (ML)] detection or interference cancellation methods 
[1-14]. In what follows, we briefly summarize the highlights of these two 
approaches as well as their advantages and disadvantages relative to the 
conventional approach wherein each user is demodulated assuming that the 
other users are not sharing the channel. A nice summary of some of these 
considerations as well as an extensive list of references on the subject can be 
found in the article by Poor [15]. 

The conventional receiver for asynchronous multiuser communication 
referred to above is comprised of a parallel bank of M (the number of users) 
receivers each of which is optimum for that particular user communicating 
over the additive white Gaussian noise (AWGN) channel alone.^ As such 
data detection for each user is accomplished independentiy of that performed 
on the other users. Although this structure is simple to implement, it is only 
optimum against the AWGN background but far from optimum against the 
highly structured nonGaussian interference produced by the simultaneous 

Us an example, if each user is transmitting an uncoded direct sequence BPSK modulation, 
then each element of the parallel receiver bank would be composed of a coherent earner 
demodulator, matched filter (chip matched filter, code despreader, and accumulator), and 

hard decision device. 
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presence of the other M - 1 users in the channel. Other disadvantages of the 
conventional scheme are: (1) it can become multiple access interference 
limited when the signal-to-background (Gaussian) noise is sufficiently large, 
and (2) it suffers from the near-far problem, i.e., high power users destroying 
the communication of low power users. The reasoning here is that when the 
received powers of the users are dissimilar, then in addition to the desired 
component, the output of each matched filter contains a component which is 
linear in the amplitude of each of the M- 1 interfering users. Thus, high 
power users contribute a disproportionate amount of interference to low 
power users making the conventional detector unable to recover their 
messages. 

In view of the above disadvantages, it became natural to ask whether 
significant performance improvement could be achieved by investigating 
receivers structures that were optimum against the background Gaussian 
noise as well as the multiuser interference. Of the multitude of contributions 
that exist in this area, the work of Verdu [1,2] is perhaps the most dted in the 
literature and the one upon which much of the other work is based. In 
Verdu's work, the receiver structme is derived based on minimizing the 
squared Euclidean distance between the received signal and the sum of the 
M asynchronous user signals, i.e, the total transmitted signal. As such, the 
presence of all M users simultaneously sharing the channel is accounted for 
in arriving at the ML receiver. The primary difierence between the structure 
that evolves from such an approach and the conventional structure is that 
joint sequence decisions are made on. the set of M matched filter outputs as 
opposed to individual bit-by-bit decisions on each matched filter output 
alone. 

While indeed such optimum multiuser algorithms offer significantly 
improved performance by alleviating the disadvantages associated with the 
conventional scheme, they unfortunately suffer from the fact that their 
complexity grows exponentially with the number of users and the length of 
the sequence. This follows directly from the fact that the optimum ML 
decision algorithm can be implemented as a d3mamic program with time 
complexity per binary decision that is 0(2^) [16]. While in many practical 
applications such performance complexity prohibits implementation of the 
Verdu algorithm, its performance is still very much of interest since it serves 
as a benchmark against which to compare other schemes with less 
implementation complexity such as those that employ interference 
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cancellation to be discussed shortly. Another disadvantage of the Verdu 
algorithm is the necessity of knowing the relative amplitudes of the various 
user signals present at the input to the receiver. One possibility around this 
disadvantage is to perform multiuser amplitude estimation [17]. An 
alternative scheme is to employ power control at the receiver which is a 
common technique used in cellular radio systems to solve the near-far 
problem. In this case, all received users are assumed to have the same power. 

The most obvious solution to the multiuser interference problem would 
be to design the user codes to have more stringent crosscorrelation properties 
since indeed if the signals were truly orthogonal this interference would not 
exist. Unfortunately, however, the near-far problem mentioned above still 
exists even for well-designed quasiorthogonal signal constellations. Thus, the 
multiuser interference problem must be dealt with and tackled from another 
viewpoint. 

One popular approach is to employ interference cancellation (IC), i.e., to 
attempt removal of the multiuser interference from each user's received 
signal before making data decisions. In principle, the IC schemes considered 
in the literature fall into two categories, namely, serial (successive) and parallel 
cancellation. With regard to the former, Viterbi [6] (see also Dent [7]) 
suggested coordinated processing of the received signal with a successive 
cancellation scheme in which the interference caused by the remaining users 
is removed from each user in succession. One disadvantage of this scheme is 
the fact that a specific geometric power distribution must be assigned to the 
users in order that each see the same signal power to background plus 
interference noise ratio. This comes about because of the fact that with 
successive cancellation the first user to be processed sees all the interference 
from the remaining M-1 users whereas each user downstream sees less and 
less interference as the cancellation progresses. Another disadvantage of this 
scheme has to do with the required delay necessary to fiilly accomplish the 
interference cancellation for all users in the system. Since the interference 
cancellation proceeds serially, a delay on the order of M bit times is required 
to complete the job. Nevertheless, Viterbi showed that the successive IC 
scheme could approach channel capacity for the aggregate Gaussian noise 
channel. As such, the scheme does not become multiuser interference 

limited. 

Parallel processing of multiuser interference simultaneously removes from 
each user the total interference produced by the remaining users accessing 
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the channel. In this way, each user in the system receives equal treatment 
insofar as the attempt is made to completely cancel his or her multiple user 
interference. As compared with the serial processing scheme, since the 
interference cancellation is performed in parallel for all users, the delay 
required to complete the operation is only a single bit time. The early papers 
that dealt with parallel interference cancellation recognized the desire to 
arrive at a structure that could be motivated by the ML approach. In 
particular, a multistage iterative approach was suggested by Varanasi and 
Aazhang [8,9] which at a given stage estimated a given user's bit xmder the 
assumption that the exact knowledge of the other users' bits in the same 
transmission interval needed to compute the multiuser interference was 
replaced by estimates of these bits from the previous stage. It was indeed this 
basic idea which led to the multistage iterative schemes subsequently 
proposed by Yoon, Kohno, and Imai [11-13] and Kawabe et al [14]. What was 
common to all of these schemes was the fact that at each stage of the iteration, 
a brute force attempt was made for each user to completely cancel the 
interference caused by all the other users. As we shall see in this paper, this is 
not necessarily the best philosophy. Rather, when the interference estimate is 
poor (as in the early stages of interference cancellation), it is preferable not to 
cancel the entire amount of multiuser interference. As the IC operation 
progresses, the estimates of the multiuser interference improve and thus in 
the later stages of the iterative scheme, it becomes desirable to cancel more of 
the interference. The motivation behind this approach can also be derived 
from maximum-likelihood considerations as was done for the brute force 
approach previously considered. 

With the above discussion in mind, this report presents a general parallel 
interference cancellation scheme that significantly reduces the degrading 
effect of multiuser interference but with a lesser implementation complexity 
than the maximum-likelihood technique of Verdu and with improved 
performance over the previously considered parallel and serial processing 
techniques. When compared with classical CDMA without interference 
cancellation, the improvement in performance is even more dramatic. Our 
scheme is suitable to the case of a nonimiform power distribution as well as a 
uniform power distribution among the users. However, we shall focus more 
on the latter since in parallel interference cancellation schemes, because each 
user sees the same amount of interference, a uniform power distribution 
achieves the best overall level of performance. In addition, although our 
scheme is suitable to asynchronous transmission, we shall assume here that 
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all users have synchronous data streams. This case results in worst case 
performance, i.e., if the data transition instants of the various users are not 
aligned, then on the average they have less of an interfering effect on one 
another. 

2.0 Multiuser Communication System Model 

We consider a CDMA communication system in which M users are 
communicating simultaneously at the same rate over a common AWGN 
channel each with a BPSK data modulation and their own pseudonoise (PN) 
code. As such the received signal is the sum of M direct sequence BPSK 
signals each with power Sf , bit time T^, and PN chip time Tq, and additive 
white Gaussian noise with single-sided power spectral density (PSD) 
w/Hz. At baseband, this signal can be written in the complex form^ 

M M 

m = ^5.(0 + n{t) = + n(t) (1) 

1=1 i=l 

M 

where for the user PN-(t) is the PN code waveform, m-(0= 

is the data modulation with bit taking on equiprobable values +1 and 
unit power rectangular pulse shape p(t) of duration T[>, and 0/ is the carrier 
phase. For the equal user power case, one would have, 5- = 5; i = 1,2,...,A/. 

We shall assume for the purpose of analysis and simulation that the users 
have purely random PN codes^ assigned to them. It is to be emphasized, 
however, that the IC schemes to be discussed in what follows apply equally 
well to any appropriate set of PN codes chosen for the users provided that the 
crosscorrelation matrix of these codes is known. In view of our assumption, 
over the zero*^ bit interval, the user's PN waveform can be expressed in 
the form 


PNiit) = '^c^p{t-kT^) ( 2 ) 

where p(0is again a unit power rectangular pulse shap>e now of duration T^, 
77 = is the number of PN code chips per data bit, i.e., the spreading ratio , 
and {cjt} is a random binary (+1) sequence. The user codes are thus specified 

^For convenience, we shall use complex notation to represent the various signals in the 
receiver. 

^For very long linear feedback shift registers, PN codes can be assumed to be purely random. 
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in terms of their normalized crosscorrelation matrix r=[y^.] where 


n = PN,mPN,(t)dr. u = 1.2.... ,M 

with 7a = 1; i = l,2,...,M. 


(3) 


3.0 Brute Force Parallel Interference Cancellation 

3.1 Tentative Hard Decisions - Equal Power, Synchronous Users 

In Appendix B we discuss several decision metrics for parallel interference 
cancellation each of which is motivated by ML considerations. All of these 
metrics suggest IC schemes in the form of iterative structures in which each 
successive stage simultaneously processes (removes multiuser interference 
from) all of the users' received signals so as to further refine the quality of 
their data estimates. The simplest of these schemes is the so-called brute 
force approach discussed in Section II of Appendix B in which an attempt is 
made at each IC stage to completely remove from each user's decision variable 
the multiuser interference contributed by the remaining users. This scheme 
is typical of what is found in the literature on parallel IC schemes [11-14]. 

We begin our discussion by considering the performance of a single stage 
brute force parallel IC scheme such as that illustrated in Fig. 1 with L = 1 
together with Fig. 2a. We further ^sume that the tentative decision devices 
associated with each user are one bit quantizers (hard decisions) as shown in 
Fig. 3a. This particular case directly corresponds to the scheme proposed in 
[11-13]. Although the results derived in Appendix B allow for arbitrary user 
powers, we assume here that all users have the same power; thus, it is 
sufficient to characterize only the performance of any one user, say the first, 
he or she being typical of all the others. Furthermore, as previously stated, we 
assume that all users have synchronous data streams and purely random PN 
codes. 


The first step is to process the received signal with a matdied filter which 
consists of despreading and demodulating^ r(f) with user I's PN code and 
carrier reference signal (both of these operations are assumed to be ideal) and 
then passing the result through a normalized I&D circuit, the output of 

^ince we are working with a baseband model, the term "remodulation" or "demodulation" 
refers to complex multiplication by the particular user's carrier phase or its complex conjugate, 
respectively. 
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which is given by 


where £* =5T, denotes 5ie bit energy, a,o is the polarity of user i's bit in the 
interval O^f^T;, n,=-L£ n(t)PN,{t)dt is a zero mean complex Gaussian 

random variable with variance £{hr} = ^o representing the thermal noise, 
= / = 2,3,...Af are the interference noises 

contributed by the other M-1 users which are modeled as independent zero 
mean Gaussian random variables each with variance 5T,.5 Also, the first 
subscript on x denotes the stage at which we are observing the I&D ou^ut 
while the second subscript denotes the particular user. This notation wiU be 
useful later on in our discussion of multiple stage cancellation schemes. The 
foregoing model of user interference as additive Gaussian noise follows from 
the assumptions made in similar analyses of CDMA systems [18,19], namely, a 
large spreading ratio tj = T, /T„ and purely random PN codes. 


Tentative hard decisions are made on the signals Zo,-; / — l,2,...,Af and are 
used in an attempt to completely cancel the other user interference. If a 
correct tentative decision is made on a particular other user’s bit, then the 
interference from that user can be completely cancelled. On the other hand, if 
an incorrect tentative decision is made, then the interference from that user 
will be enhanced rather than cancelled. A quantitative description of this 
will be given when we model the signal upon which final decisions are 
made. As we shall see, the performance analysis associated with this model is 
complicated by the fact that the tentative decisions are not independent of one 
another. More about this shortly. 


Defining the Gaussian noise variable and PN crosscorrelation for the 
ith user analogous to the above definitions for user 1 (see (A.4) and (A.6) of 
Appendix A), then after respreading/remodulation, interference cancellation, 
and despreading/demodulation, the normalized output of the I&D 
corresponding to the final decisions is given by 


^The normalized interference noises i — 2,3,.. .,Af have variance equal to the 
reciprocal of the spreading ratio, i.e., TJ 
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•*u 



i>2 


where 


(5) 


A = flio - sgn Rej VIJ fl,o + 


is a three-valued (0, +.2) indicator random variable whose magnitude 
represents whether or not a correct tentative decision is made on the 
user's bit. It is tempting to model the ^'s as independent random variables. 
Unfortunately, this leads to optimistic results (when compared with the true 
performance results obtained from simulation). In addition to the fact that 
the ^,'s are not themselves independent, they are also dependent on the PN 
crosscorrelations, i.e., the y„’s. Fortimately, however, the )3/s are not 
strongly dependent, i.e., aside from the weak dependence between the 
Gaussian noises n, and rij (see (A.5) of Appendix A), the only signal terms 
that preclude complete independence of say and j8y are in )3/ and 
^ioYji=^io7^ in A/- Hence, for sufficiently large M, it is reasonable to assume 
a Gaussian model for the total residual (after cancellation) interference term 
/j in (5). The accuracy of this model will improve as M increases (actually as 
the munber of nonzero terms in 7j increases which implies a high tentative 
decision error rate). Later on we shall numerically compare the performance 
results derived from this analytical model with those obtained from a true 
computer simulation of the receiver to establish the accuracy of this model. 


Assuming then a Gaussian model for /, (note that /, is not zero mean), 
then the average probability of error associated with the final decisions is 
given by 


n(« = ipr{Rc{x„ >0|,_..,}}+ipr{Re{x„ <o|,..,}} 


= > o|,..-i}} = rr|A', > -*) 


where^ 


simplify the notation here and in what follows, it is understood that the statistical 
mean cos((p- — ) is computed under the hypothesis = — 1. 
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= Re{n,e"'^ + /, -A} = + ^ITX AriiCOs(^i -^)- V^TX Ayi.cos(^. - ^) 

i -2 «-2 

( 8 ) 

is the effective noise seen by user 1 after cancellation which in view of the 
above is modeled as a real zero mean Gaiissian noise random variable whose 
thermal noise component has variance o^=N^ /2. It is straightforward to 
compute the variance of as 

<tI=N,/ 2 + E,(M - DPfrl cos^(^. - ^0 - £*(Af - l)*(Aru cos(^. 

+E, (M - l)(Af - 2)p,YuPjrij cos(0i - ^ )cos(^; - ^ ) (9) 

+2V^(A/-l)iV,Ay.,cos(0, -^) 

where i can take on any value from the set 2, 3,..., M. Hence, from (7), the 
average probability of error can be obtained as 


P,(E) = 



where 






l + 2^(M-l) +(A/-2)44,+^iV^ 


( 10 ) 

^i.=Ari;COs(0i-^,) 

( 11 ) 


is an SNR degradation factor (relative to the performance of a single BPSK 
user transmitting alone) and Qix) is the Gaussian probability integral 
defined by 


2}^ 

Thus, the evaluation of P*(£) reduces to the evaluation of the various 
statistical averages involving required in (11). These statistical averages, 
which must be performed over the Gaussian noise and interference random 
variables as well as the uniformly distributed carrier phases, are not trivial to 
compute in closed form. The primary difficulty lies in the averages that must 
be performed over the interference crosscorrelations, y^, which from their 
definition in (3) are binomially distributed in the interval (-1,1). In particular. 
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for the assumption of random user codes, 
that takes on values \-2klr], ife = 0,1,2,.. .,r/ 


is a discrete random variable 


with probabilities P* = 




further complication lies in the fact that although the y^'s are pairwise 
uncorrelated (see Appendix E), they are not independent. All is not lost 
however. Since the variance of y^ is equal to \IN^ = llr\, then for large t] 
(the typical case of interest), the p.d.f. of y^ will be quite narrow with respect 
to its domain and hence for analysis purposes we can approximate the y,^'s as 
being Gaussian distributed. Also in view of this Gaussian assumption on the 
Yy’s, we can now further invoke their independence becavise of their pairwise 
uncorrelated property. Even with the assumption of independent Gaussian 
random variables for the y^'s, evaluation of the statistical averages required 
in (11) is still quite complex. Nevertheless, with the help of some useful 
closed form integrals (see Appendix F), the needed moments can be obtained 
in the form of definite integrals of tabulated functions with the following 
results:^ 


r=-J_r 


2cr^ cos^^ 

aVcos^^ \ 




expe- 


ct 


2(l + aVcos^^) 


\d<p (13a) 


^ 1 [T 4 4 / 

q = — .—a cos A 

2;rJ« ' 


1 f2* 


a‘ 


x3/2 


1 + aVcos^^ 


cxp^ 


a 


2(l + aV*cos*0)J 


d(P 


(13b) 




^Despite the Gaussian assumption on y, the integration on this random variable in (13c) 
cannot extend from — ®o to «» . The reason for this is that the normalized correlation, p, 
which is directly proportional to y in accordance with (14) must be restricted to the 
interval (-1,1). Hence, in the evaluation of (13c) we actually replaced the — and «> limits 
for the integration on yby the limits -1 and 1 thereby treating it as a truncated Gaussian 
random variable. 
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where denotes the required bit energy-to-noise spectral density ratio 

for M users connmunicating simultaneously, each of which operates at an 
average bit error rate P^iE). 


It is common in analyses of CDMA systems [18] to define a degradation 
factor, D as the ratio (in dB) of the EJN^ required to achieve a given bit 
error rate in the presence of M users, namely, (E^/No), to that which would 
be required to achieve the same level of performance if only a single user was 
communicating, namely, {EJNo\. By the definition of (E^/Z/o),/ l^^ve 
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To obtain the degradation factor for a given value of we substitute 

Dx(£./W.),=Z)x[i[Q-(/l(£))r] for (E./JV.) ^ in (14) which in turn is 

substituted in (13). Then using the given value of /*(E) one can solve for D. 
Unfortunately, a closed form expression for D cannot be obtained so the 
results will be obtained numerically. Before presenting these numerical 
results, however, we briefly review the analogous results for conventional 
CDMA and the successive (serial) interference cancellation sdieme proposed 
by Viterbi [6] (later patented by Dent [7]) since we shall use these as a basis of 
comparison to demonstrate the increased effectiveness of parallel 
cancellation. 


3.1.1 Comparison with Conventional CDMA and Successive Interference 
Cancellation 

In a conventional CDMA system, there is no attempt made to cancel the 
other user interference. Hence, (£’ 4 /A/<,)j is given by 

\ §b = 

Thus, the degradation factor, D, is [18] 


P _ 1 


(17) 


For the successive cancellation scheme [6], Viterbi showed that to 
guarantee that each user in the system sees the same amount of interference 
from the other users, the user powers should be assigned as 


5 * = 5 . 


Nn 




jk = Af,Af-l,...,2 


(18) 


where is the power of the user to be processed last (the weakest one) and S^f 
is the power of the user to be processed first (the strongest one). Distributing 
the powers eis in (18) ideally guarantees that all users see the same ratio of 
signal power to effective noise spectral density and thus the user to be 
processed first (the one that sees all the user interference) is not any SNR 
disadvantage relative to the user to be processed last (the one for which all 
interference has been removed). In view of the above, the degradation factor 
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for the fcth user is given by 


( 19 ) 


where {E^/N^)j^ denotes the required bit energy-to-noise spectral density ratio 


for the user. The average degradation factor, D for the M user system is 
obtained by averaging (19) over k which yields 


If.^ (l + n-'(£./A^.),)“-l (20) 

uh ‘ Mn-'(EJN,\ 

It should be emphasized that the result in (20) ignores the effect of decision 
errors made at the various successive interference cancellation stages, that is, 
the interference cancellation is assumed to take perfectly. As a result, 
numerical results derived from (20) will be optimistic when compared to the 
actual performance of the scheme. 


3.1.2 Numerical Results 

To illustrate the significant performance advantage of the brute force 
parallel interference cancellation scheme in Fig. 1 combined with Fig. 2a, we 
consider a plot of D versus M for an average bit error probability® 
/>^(£) = 10"^ and a spreading ratio rj = l00. Fig. 4 shows the analytical 
performance of the three schemes (conventional, successive interference 
cancellation, parallel interference, cancellation) as well as computer 
simulation results for the latter. The analytical models that describe the 
simulation are presented in Appendices A and B. We see that for the 
conventional and parallel interference cancellation schemes there exists a 
user capacity limit in that regardless of how much one is willing to increase 
(£fc/No)j, (for a given (£*/A^o),/ or equivalentiy, a given P^(E)), the required 
bit error rate cannot be achieved if more th<m Af^ users simultaneously 
access the system. For conventional CDMA 


=1 + 




A 




( 21 ) 


whereas for the parallel interference cancellation scheme the solution is 

®The value of /**(£) = 10"^ is chosen to allow for obtaining computer simulation results in 
a reasonable amount of time. 
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determined from 


10-" = Q 


l + 

- 1)[S - - l)(^f + < - 2)S] ^ 

together with the moments in (13) where now 


( 22 ) 


It is emphasized that the user capacity limit for the parallel interference 
cancellation scheme comes about entirely because of the finite probability of 
error associated with the tentative decisions. From Hgure 4 it appears that 
the successive interference cancellation does not have a user capacity limit. 
This is because in [6], it was assumed for this scheme that the interference 
cancellation is perfect, i.e., the effect of decision errors at the various 
interference cancellation stages were not accounted for. 


Comparing the analytical and simulation results for the parallel 
interference cancellation scheme, we observe that the two are virtually in 
perfect agreement. Any discrepancy between the two stems from the analysis 
assumptions made in connection with (5) of: a) an analytical Gaussian model 
for the total residual user interference, and b) independent Gaussian code 
crosscorrelations (y,y's), whereas the computer simulation makes no such 
assumptions and thus predicts the exact performance. For smaller user error 
probabilities, e.g., Pj(£) = 10"*, one is likely to expect a larger discrepancy 
between the analytical and simulation results since, in this case, the number 
of nonzero (residual) interference terms could be small in comparison with 
the number of users and hence the Gaussian residual interference 
assumption might well become weak. Under these circumstances, one would 
be best to rely entirely on the simulation results. 

3.2 Tentative Hard Decisions - Unequal Power, Synchronous Users 


The results of the previous section can be generalized to the case where the 
users have imequal powers, i.e., = Let a^ = SJSj denote the ratio 

of the power of the ith user to that of the ;th user who is arbitrarily 
considered to be the desired user. After interference cancellation, the 
normalized output of the I&D corresponding to the final decisions of user j 
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is by analogy with (5) 


«-l 

i»j 

(24) 


where n- =-^ r*n(0/*A^, (0^.y = 1.2 M is a zero mean complex Gaussian 

random variable with variance £|jny[*|=A^o representing the thermal noise of 


the /th user, y^-=— J *£iV.(t)FlV,(f)dr, / y is the normalized interference 

produced by user i on user / [see (3)], and £^^=5/* is the bit energy of the 
user. Also, analogous to (6), is now defined by 




f 

= fl;o-sgn 

Re-i 
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msl 


•n.e 




(25) 


Following steps analogous to (7) - (9) we arrive at the desired result for the bit 
error probability of the desired (the ;th) user, namely. 


where 


N ‘ 

\\ ‘'^0 j 


(26) 
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^yj=Ay>cos(^i-^;) 

( 27 ) 

Clearly for all a. 's = 1 and; = 1, (27) reduces to (11). 
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As an example, consider a group of M users with powers exponentially 
distributed (linearly distributed on a dB scale) over a range of 10 dB between 
the minimum and the maximum. This model might correspond to a 
distribution of mobile users that are exponentially distant from the base 
station within a cell. Assume that we fix the error probability of the lowest 


power user (assumed to be user 1 for convenience of notation) equal to lO"*^ 
(obviously, all others would then have lower error probability.) Then, Fig. 5 
illustrates the degradation factor, D^, of user 1 versus M. For comparison, 
the results corresponding to conventional CDMA with the same user power 
distribution are also shown in this figure. These conventional CDMA results 
are obtained from Eq. (9) of [18] which is a generalization of (17) for the 


unequal power 


case, namely, Dj = 





By comparison 


with Fig. 4, we observe that in the imequal power case, parallel interference 
cancellation offers more of a relative advantage compared to conventional 
CDMA than in the equal power case. The reason behind this observation is 
that the larger power of the other users (which are producing the multiuser 
interference to user 1) produces tentative decisions with a smaller error 
probability which in turn results in a better degree of cancellation with regard 


to the final decisions. 


3.3 Null Zone Tentative Decisions - S)mchronous Users 

Much like the idea of including erasures in conventional data detection to 
eliminate the need for making decisions when the SNR is low, one can 
employ a null zone hard decision device (see Fig. 3d) for the tentative 
decisions in Fig. 2a to further improve the fidelity of the interference 
cancellation process. The idea here is that when a given user’s signal to 
interference ratio is low, it is better not to attempt to cancel the interference 
from that user than to erroneously detect his data bit and thus enhance his 
interference. Following the development in Section 3.1 for a single stage 
scheme with equal power s)mchronous users, then the normalized output of 
the I&D corresponding to the final decision on user I's bit o,o is still given by 
(5) with Pi now defined by 




Re-< 



V 








m-l 

m#i J 


A=^.o-nsgn 


(28) 
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where "nsgn" denotes the null zone signum function defined by 


( 29 ) 


X x>C 

nsgn z = '0, 

-I, x<-C 

Here Pi takes on possible values (0, ±1, +2) and its magnitude is an indicator 
of whether a correct decision is made (ith user s interference is perfectly 
cancelled), no decision is made (ith user's interference is unaltered), or an 
incorrect decision is made (ith user's interference is enhanced). Once again 
making a Gaussian assumption on the total residual inteference, then since 
the final decisions are still made as hard decisions, the average bit error 
probability is still given by (10) together with (11) with now 
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where 


ir±fD o A A X /» A +^2 ±pfiig2g^(cos"^ +cos"^) 

^ ^ l-p^a^B^B2Cos^p^cos^^2 

G(B,.fi2,0p02)= J _ p2^4^ cj2^^ C052^^ 


(30e) 


with fli,52 defined as before in (14) and C = C/V^ ^ normalized decision 
threshold which should be chosen to minimize D for a given Pj(£) and 
(£^/7Vo), determined from (15). 


Numericad results for the case where a null zone limiter is used as the 
tentative decision device are illustrated in Fig. 6 together with the 
corresponding results previously given in Fig. 4 for the case of hard limited 
tentative decisions. For the specified processing gain and average bit error 
probability (the same values as in Fig. 4), we see that using a null zone limiter 


19 


allows the maximum number of users that can be supported to be increased 
by about 10% - a modest improvement in performance. The normalized 
threshold has been fixed at ^' = 0.4 which appears to be the value that gives 
the best performance improvement. For an unequal (exponentially 
distributed) power distribution among the users, the corresponding results 
using null zone tentative decisions are superimposed on the previous 
discussed results illustrated in Fig. 5. Here, the nonnaUzed threshold^ has 
been fixed for each user at ^ = 0.25 which appears to give the best overall 
performance improvement for the specified user power distribution. Again 
we see a modest improvement in performance. 

3.4 Multiple Stage Interference Cancellation 

The single stage scheme of Fig. 1 combined with Fig. 2a can be improved 
upon by cascading multiple stages (L > 1 in Fig. 1) of parallel interference 
cancellation. ^0 The idea here is to repeatedly improve the fidelity of the M 
tentative decisions since each successive stage sees less and less interference. 
Note that in principle this idea is similar to what Viterbi accomplishes in the 
serial interference cancellation scheme except that here at each stage we 
simultaneously act on the interference from all M users rather than one user 
at a time. An analysis of the performance of such a multistage scheme is 
difficult if not impossible to obtain due to the fact that the tentative decisions 
at the interference cancellation stage depend on the tentative decisions at 
the (j-l)®* stage. Because of this difficulty, numerical results for the 
performance of the multi-stage parallel interference scheme must be obtained 
from computer simulation. Illustrated in Fig. 7 are performance results for 
two and three stage parallel IC schemes with hard tentative decisions, 
respectively. The corresponding one stage results are repeated on these 
figures for the purpose of comparison. We observe that there is significant 
gain to be achieved by going to more than one stage. 


^In the unequal user power case, the normalized threshold for the ith user is defined by 
C = C. / where C, is tbe actual decision threshold for this user. Thus, a fixed 
normalized threshold implies a different actual threshold for each user in accordance with 
his or her energy (power). 

For the case where the tentative decisions associated with each user are one bit 
quantizers (hard decisions) the scheme again reduces to that proposed in 111-13]. 
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3.5 Linear (Infinitely Quantized Soft) Tentative Decisions - Equal Power, 
Synchronous Users 

One disadvantage of the parallel cancellation scheme with hard or null 
zone limiter tentative decisions is that, in order to jjerform the respreading 
and remodulating operations, the receiver must ideally have complete 
knowledge of each user's power, carrier phase and frequency, and PN code 
chip timing epoch. Since in practice the receiver does not have knowledge of 
these parameters, it must estimate them. One simple way of drciunventing 
some of these problems is to use linear (infinitely soft quantization) tentative 
decisions (see Fig. 2c). Since the signal component of the output of the 
tentative decision devices is now linearly proportional to the square root of 
the user powers, it is no longer necessary for the receiver to estimate these 
powers prior to the cancellation operation and thus the ‘JeH gains following 
these devices in Fig. 2a may be eliminated. Another simplification that is 
now possible is that the final decisions can be performed with a differential 
(rather than coherent) detector thus eliminating the need for carrier 
synchronization at all stages. There is of course a performance penalty 
associated with using differential rather than coherent detection; however, it 
was shown in [20] that by using a multiple symbol (as opposed to the 
conventional two-symbol) differential detector one can fill the gap between 
these two performances. The primary disadvantage to using linear tentative 
decisions is that additive thermal noise from each user is now directly 
introduced into the interference cancellation process. This will result in a 
performance that is inferior to the .hard and null zone tentative decision 
schemes but still better than conventional CDMA which employs no 
interference cancellation at all. Furthermore, one can now in principle 
analytically compute the performance of a multiple stage parallel interference 
scheme although the analysis becomes quickly complex as the number of 
stages increases beyond two or three. 


We first consider the analysis of a single stage scheme using linear 
tentative decision devices and with equal power, synchronous users. 
Analogous to (4), the normalized output of the ;th I&D circuit in Fig. 2a is 

given by 




( 31 ) 
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After respreading/remodulation, interference cancellation, and 
despreading/ demodulation, the normalized output of the I&D corresponding 
to final decisions of user 1 (again assumed without loss of generality to be 
the desired user) is given by 


•^1 ““ ^0 


M 
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j>2 


i«2 


km2 

k^i 


M 

I 

i-2 




(32) 


We note that the decision statistic for user 1 now includes an attentuation of 
the desired signal component as well as additional thermal noise due to the 
presence of the other M-1 users. The middle term in (32) represents the 
residual user interference after cancellation. 


Exact evaluation of average error probability would be obtained by 
computing the bit error probability conditioned on a given set of user codes 
(random sequences) which takes the form of a Gaussian Q-function and then 
averaging the result over all user codes. An approximation to this is obtained 
by treating the interference as an additional Gaussian noise source thus 
allowing a solution for the average error probability in the form of a single 
Gaussian Q-function. This approach is, in principle, the same as that 
previously used for the nonlinear IC case except that here we can exactly 
evaluate the first and second moments of the interference using the true 
binomial distribution of the y/s. In particular, the mean (conditioned on 
Ojp = 1) and variance of x,, are computed in Appendix D. Using these results 

the average bit error probability for hard decisions made on =Re{x,,} is 


P,(E) = 



(33) 


where cr,^, denotes the variance of Xg^^ and 
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or, equivalently, the degradation factor is given by 

M-l (A/-l)(Af-2) 

V rf (35) 

(f * /^o). ^ Y JeA( + 2 + 

I V ) UoJA A V } 

Analogous to (21), the maximum number of users, that can be 

supported with a single stage linear cancellation scheme is obtained by setting 
D = oo oT equivalently setting the denominator of (35) equal to zero. An 
explicit solution for is cumbersome. For large M, one has the 

asymptotic (as becomes sufficiently large) solution 



Thus, ignoring the "1+" term on the right hand side of (21), (valid for large 
A/„,), we see that the maximum number of users that can be supported by a 
CDMA system using a single stage of linear interference cancellation is 
increased relative to that of a conventional CDMA system with no 
interference cancellation by the factor in the limit of large {E^/Nq)^. 

Figure 8a is an illustration of D versus M obtained from (35) for rj = 1(X) and 
/^(£) = 10“^, i.e., the same processing gain and average bit error probability as 
in Fig. 4. Superimposed on this figure are the results corresponding to 
conventional CDMA obtained from Eq. (17). We see that for low values of 
D, improvement is obtained with the brute force linear cancellation method 
relative to conventional CDMA whereas asymptotically for large D (i.e., 
where multiuser interference dominates), we indeed obtain a slight loss in 
performance. Clearly then (£^7^^)^ corresponding to /^(£) = 10“^ is not 
sufficiently large to justify the use of (36). Furthermore, the brute force linear 
cancellation method is quite inferior to the brute force hard-limited tentative 
decision case (see Fig. 4). In Figure 8b are the analogous results for 77 = 100 
and a smaller error probability, namely, i^(£) = 10~*. Superimposed on this 
figure are again the results corresponding to conventional CDMA. We now 
see that an improvement relative to conventional CDMA exists for both 
small and large values of D. Still, however, the linear cancellation method 
is far inferior to the nonlinear cancellation method using hard tentative 
decisions. 
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Before proceeding to a discussion of multiple stage linear IC schemes, we 
wish to compare the "exact"^^ degradation factor found in (35) with that 
which would be obtained by invoking an independent Gaussian assumption 
on the y^'s as discussed in the latter part of Appendix D. Using (D.3‘) - (D.7) 
in place of (D.3) - (D.7), the degradation factor becomes^2 

^ M-1 

V (35') 

{E,/N,\ f f£.'|r (M-lXA^ + 2)' | 

I V J W)i J 

Comparing (35’) with (35) we see that the results are similar except for higher 
order terms in rj-\ To quantify this comparison, the degradation factor, D, 
as computed from (35') is plotted in Fig. 8 along side the "exact" result 
computed from (35). Also shown here are the simulation results which 
correspond to the true exact performance. We observe that the degradation 
factor computation based on the independent Gaussian assumption for the 
y^'s is reasonably close (but somewhat optimistic) to the "exact" result based 
on binomially distributed y./s that are not necessarily independent. The 
reasonable similarity of the "exact" and approximate degradation factor 
results for the single stage linear IC scheme will serve as a justification for 
using the latter in making comparisons among the performances of multiple 
stage linear IC schemes where exact computation of D is extremely tedious. 

As discussed in Section 3.4 matters can be improved by going to a multiple 
stage interference cancellation configuration. For two stages of cancellation 
with linear tentative decisions as illustrated in Fig. 1 combined with Fig. 2a, it 
can be shown that the normalized output of the I&D corresponding to the 
final decisions of user 1 is (by analogy with (32)) 


”ln referring to the degradaHon factor in (35) as "exact" we mean that the true binomial 
(as opposed to a Gaussian) distribution of the y/s is used in computing the moments of these 
random variables. We hasten to point out, however, that the degradation factor in (35) is 
still not truly exact (i.e., it will not totally agree with the simulation results) because, as 
mentioned in the text, the total residual interference is still being analytically modeled as a 
Gaussian random variable [see the discussion preceding Eq. (33).] 
l^Note that (35 ) results in the same asymptotic as (35) does, namely, (36). 
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The first term again represents the desired signal component, the middle 
term the other user interference and the third term the additive thermal 
noise. Using Eq. (E.3) of Appendix E, it is straightforward to show that the 
mean of the first term (conditioned on flio = D is equal to 
.y^l + (Af-l)(Af-2)/77^]while the other two terms ii;i (37) have zero mean. 
Thus, for the two stage linear cancellation scheme 
£{;c 2 i|^^.,} = V^l + (^“i)(^”2)/77*] which represents a signal enhancement 
caused by the presence of the other users. The variance of each term can, in 
principle, be computed analogous to Eqs. (D.3), (D.4), and (D.5) of Appendix D 
but with considerably more difficulty which makes this approach impractical. 
Thus, we shall proceed as in the latter part of Appendix D by invoking an 
independent Gaussian assumption for the y/s. When this is done, the 
following results are obtained. 


The mean of the desired signal term (conditioned on = 1) is now equal 
to -Je^. Similarly, the means of the other two terms in (37) are still equal to 
zero. The variance of the self noise of the desired signal component, user 
interference, and equivalent thermal noise^^ are given by 
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(38) 


(39) 


^^Recall that the noise components {nj are conditionally correlated as per Eq. (A.5) of 
Appendix A. 
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FinaUy, combining (38H40) we have analogous to (34) 
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or solving for the degradation factor 
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Analogous to (36), the maximum number of users that can be supported with 
a two stage linear cancellation scheme is obtained by setting D = oo in (42). 
which for large can be approximated by 



V 1^0 

Thus, again ignoring the ’T+" term on the right hand side of (21), the 
maximum number of users that can be supported by a CDMA system using a 
two stage linear interference cancellation is increased relative to that of a 
conventional CDMA system with no interference cancellation by the factor 

^[(£j/No)if • An illustration of D versus M obtained from (42) is super- 
imposed on the results of Figures 8a and b. We observe that the addition of a 
second stage of cancellation gives us a modest improvement in performance. 


Obtaining performance results for more than two stages of linear 
interference cancellation can in principle be accomplished but requires 
formidable analytical evaluations. The procedure for K stages and the 
difficulty in obtaining closed form results are described as follows. Let 
^ -{y r I 1 = 012 K denote the column vector of M I&D outputs 
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at the ith stage of the interference canceller. Furthermore, let A = F-I where 
r = | 7 ^.| is the M'x.M matrix of normalized PN crosscorrelations and I is the 
M X M identity matrix. Since = 1 for all i = and = y^, then A is 

a symmetric matrix with zeros along its diagonal. In tenns of these 
quantities, the output of the i th stage is given by the iterative relation 

X. = Xo - d>*AOx,._, (44) 

where 

* 0 = 

a = (fl,o,fl 20 — and n = (/ii,/j 2 ,...,nj|,)’^ respectively denote the data and 
noise column vectors, and ^ is the complex carrier phase diagonal matrix 
defined in Eq. (A.12) of Appendix A. The solution to (44) for the output of the 
fCth stage can be easily shown to be 

Xjt = - •>/^0*(-A)*^‘"'4>a + <I»*^(-A)‘n (46) 

j>0 

Defining the column vector u = (1,0,0,. ..,0)^, then the output of the Kth stage 
for user 1 is 

X,, = u^’xjf = - ^u^O‘(-A)^*'«>a + u^<I»*5^(-A)'n (47) 

4»0 

Letting a = OjoU + a where a = , then (47) can be rewritten as 

Xjfi =u^Xjf =.y/^Oio(l-u^<[>‘(-A)*'"4»u)-.^u’’4>*(-A)*^'<ta‘ +u^O*2^(-A)‘n 

i«i0 

(48) 

Note that for an square matrix C = |c^|, we have u^Cu = Cjj. Thus, since the 
diagonal elements of 0*(-A)*^^'C> are identical to the diagonal elements of 
(-A)*^"^’, we can simplify (48) to 

- u^(-A)'*'u)- Vj^u’’0’(-A)'*‘0a' +u’’0*]£(-A)‘n 

««0 

(49) 

For X = 1 and X = 2 , it is straightforward to check that (49) reduces to (32) 
and (37) respectively. 


The first term in (49) is real whereas the second and third terms are 
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complex for all values of K. To compute the SNR of user 1 after K stages of 
interference cancellation we need the mean and variance (conditioned on 
Oio = 1) of the real part of the decision variable in (49). These are given by 

£{Re{a:jr, = l} = - u’‘(-A)*^*u] (50) 

and 

= '} = £.{(“^(-A)"'*"u)‘ -[u^(-A)‘*'u]’ I 

= |.|u'(-Af ‘-'V +(u^-A)'‘*"uf -2[u'eATuf } 

2 >*o ' 

(51) 

where for convenience of notation the overbar denotes statistical expectation. 
Also note that 

ii(-A)‘*'= iy+;)P^+ (52) 

isO 1*0 

I even i even 

Squaring (48), dividing by (49) and equating the result to 2{E^INq)^ gives a 
relation from which the degradation factor D can be solved, namely, 

£) = ‘“O r±i — 

(l - uf - f|- j {u'(=lp^u 4- (u'(-A)'“»uf - 2[u'(-A)‘*'u]’} 

(53) 

or using (52) 

( K , 2K+\ 

X(-A)‘ - X “ 

_ l«0 jmK*\ J 

"" (l - u^(=Apu)’ - f|s- 1 + (u'(-A)“*"u)‘ - 2[u'(:7r uf} 

(54) 

Substituting as before ^[Q~'{P^(E))f for {EJN^X allows D to be evaluated for 
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a given value of /*(£). 

Although in principle either an "exact" or approximate expression for D 
can be obtained from (54) depending on whether the elements of the matrix 
r are modeled as binomial or Gaussian, the difficulty in arriving at a closed 
form result for D lies in the evaluation of the first element of the mean of 

A', i integer and also the mean of the square of the first element of A*. Using 
Mathematica software and a Gaussian model for the y^'s, we have been 
successful in carrying out these evaluations for X = 1, 2, and 3 leading to the 
following results: 


u^(-A)*u = 0, n odd 



(55) 

For X = 1 and 2, using the results of (55) in (54) gives the closed form results 
for D found in (350 and (42) respectively. For X = 3, the corresponding 
result is given by 



An illustration of D versus M as computed from (50) is superimposed on 
the resiilts of Figs. 8a and 8b. We observe fi-om Fig. 8a that for P^{E) = 10"*, the 
performance of the three stage linear cancellation scheme is actually worse 
than that of the two stage scheme. A possible explanation for this seemingly 
strange result is as follows. Observe that the first term in the denominator of 
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D in (54) represents the normalized square of the mean of the decision 
variable [see (50)]. The part of this term correponding to u^(— A)*^^*u 
represents self-noise which in view of the first equation in (55) is absent when 
K (the number of stages) is even and present when K is odd. Thus the the 
first term in the denominator of D is equal to unity when K is even and less 
than unity when K is odd, i.e., this term alternates as a function of the 
number of stages between unity and something smaller than unity. If 
{EJNo\ is small enough, as would be the case for high user error probability, 
so that the first term in the denominator of (54) tends to dominate, then for a 
given number of users the degradation will alternately get better or worse as a 
function of the number of stages. This altematation has been observed by 
simulation for 1, 2, 3, and 4 stage IC schemes, i.e., 2 stages is better (smaller D) 
than 1 stage, 3 stages is worse than 2 stages, but 4 stages is again better than 3 
stages but by a smaller amount. At small user error probabilities where 
{El,/ No)i is isrge, this alternation does not take place. To demonstrate this. 
Fig. 8b is the corresponding plot to Fig. 8a for /**(£) = 10”*. We observe here 
that the performance of the three stage linear cancellation scheme is better 
than that of the two stage scheme as one might ordinarily expect to be the 
case. Continuing beyond X = 3 is extremely tedious analytically and is not 
justified from a practical standpoint in view of the limited amount of 
performance improvement anticipated. 

4.0 Improved Parallel Interference' Cancellation Based on Joint Maximum- 
Likelihood Considerations 

In the previous sections, we discussed a so-called brute force multistage IC 
scheme in which the receiver attempts at each stage (iteration) to fully cancel 
the multiuser interference. In Appendix B it is cwgued that, since in the early 
stages of interference cancellation the tentative decisions are less reliable than 
they are in later stages, the brute force philosophy of entirely cancelling the 
interference at each iteration stage may not result in a receiver with the best 
performance. Indeed, it is suggested there that a better philosophy might be 
one which in the early stages only partially cancels the multiuser interference 
with the amount of cancellation increasing as one continues to iterate toward 
the ultimate final data decisions, i.e., as the fidelity of the tentative decisions 
improves. With this philosophy in mind, a decision statistic is derived in 
Appendix B based upon joint ML-type considerations which leads to the 
multistage receiver structure of Fig. 1 with fcth stage as in Fig. 2b. The 
modification of Fig. 2a which results in Fig. 2b is the inclusion of a parameter 
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Pt to aUow for partial cancellation of the multiuser interference at the fcth 
stage. Intuitively, one would expect that the value of p* (which depends on 
the particular stage through the subscript k) would monotonically increase as 
one progresses toward the final data decision. Indeed, the numerical results 
to be presented later on bear out this intuition. 

We begin our discussion of this modified parallel IC cancellation scheme 
by once again considering the case of tentative hard decisions with equal 
power, synchronous users. 

4.1 Tentative Hard Decisions - Equal Power, Synchronous Users 


For a single stage scheme, it is straightforward to show that the 
normalized output of the I&D corresponding to the final decisions is still 
given by (5) where now takes on the more general form 
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Note that for p, = 1, (57) reduces to (6). Following the development in Section 
3.1, it is straightforward to see that the average probabiUty is stiU given by (10) 
together with (11) where A (57) rather than (6) is used in the definition of 
When this is done, the required moments of (11) can be determined by 
applying the following replacements to the moments in (13). 


•all ^ Pi 


(58a) 



2 



(58b) 

(58c) 


2 


^ Pi 



(58d) 


which reduce, respectively, to (13a-d) when p, = 1. For Pi - 0, we have 
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= ^ = 0 and ^=i.=< 7 V 2 = l/27j. Thus, for this case, A of (11) 
simplifies to 

1 (59) 

1 + (M — 1)7] '(£j/No)j{ 

Since from (10) and (15), (f*/No)i = then from (59) and the 

definition of the degradation factor, D, we have 

D - = J_ = 1 + (M - l)T7"‘(^fc /^o)j, = 1 + (A^ - (60) 

or solving for D, 

ry I (61) 

which agrees with (17), i.e., the degradation factor corresponding to 
conventional CDMA with no interference cancellation. 

Figure 9 is a plot (obtained by computer simulation) of the degradation 
factor D versus the number of users M with p, as a parameter and the same 
values of P,{E) and t] as in Fig. 4. We observe that an optimum (in the sense 
of maximizing M for a given D) a exists whose value is other than unity. 
That is, for a spedficed amount of multiuser degradation, partial interference 
cancellation allows a considerably larger number of users to be supported 
than the number permitted by the brute force IC scheme of Section 3. 

4.2 Tentative Hard Decisions - Unequal Power, Synchronous Users 


The results of Section 4.1 can be generalized to the case where the users 
have unequal powers in the same manner that the results of Section 3.1 were 
generalized to those of Section 3.2. In particular, the normalized output of 
the I&D corresponding to the final dedsions of user j is still given by (24) 
with A I'O'v taking on the more general form in (57). Hence, the bit error 
probability of the desired (the ;th) user is stiU given by (26) together with (27) 
where, however, the moments required in (27) are obtained from (58a-d). 

4.3 Null Zone Tentative Decisions - Synchronous Users 


The results of Section 4.1 can be generalized to the case where a null zone 



32 


hard decision device (see Fig. 3b) is employed for the tentative decisions in 
Fig. 2b in the same manner that the results of Section 3.1 were generalized to 
those of Section 3.3. Following the development in Section 4.1 for a single 
stage scheme with equal power synchronous users, then the normalized 
output of the I&D corresponding to the final decision on user I's bit is still 
given by (5) with now defined analogous to (57) by 
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where "nsgn" again denotes the null zone signum function defined in (29). 
Note that for p, = 1, (62) reduces to (28). Following the development in 
Section 3.3, it is straightforward to see that the average probability is still 
given by (10) together with (11) where of (62) rather than (28) is used in the 
definition of When this is done, the required moments of (11) can be 

determined by applying the moment replacements of (58a-d) to (30a-d). 

A plot of the degradation factor D versus the niunber of users M with 
as a parameter would again reveal that an optimum (in the sense of 
maximizing M for a given D) pj exists whose value is other than unity. 
That is, for a spedficed amount of multiuser degradation, partial interference 
cancellation again allows a considerably larger number of users to be 
supported than the number permitted by the brute force IC scheme of Section 
3. 

4.4 Linear (Infinitely Quantized Soft) Tentative Decisions - Equal Power, 
Sjmchronous Users 

As in Section 4.4, we first consider the analysis of a single stage scheme 
using linear tentative decision devices and with equal power, synchronous 
users. Analogous to (32), the normalized output of the I&D corresponding to 
the final dedsions of user 1 (again assumed without loss of generality to be 
the desired user) is given by 
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(63) 

As in the brute force IC case, the dedsion statistic for user 1 includes an 
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attentuation of the desired signal component, residual user interference after 
cancellation, and additional thermal noise due to the presence of the other 
M“1 users. The mean (conditioned on flio=l) variance of jTh can be 
evaluated in a form analogous to the results given in Appendix D. Using the 
results of these evaluations, the average bit error probability for hard 
decisions made on the real part of Xu is again given by (33) where now 
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or, equivalently, the degradation factor is given by 
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(65) 


For p, =0, this result reduces to (17), i.e., the degradation factor corresponding 
to conventional CDMA with no interference cancellation. Also, for Pi = 1, 
(65) reduces to (35) corresponding to the case of full (brute force) interference 
cancellation. Figures 10a, b are plots' of D versus M obtained from (65) with 
Pi as a parameter and the same processing gain and average bit error 
probability as in Figures 9a, b. We observe that, as in the hard tentative 
decision case, an optimum (in the sense of maximizing M for a given D) Pi 
exists whose value is other than unity. That is, for a specificed amount of 
multiuser degradation, partial interference cancellation again allows a 
considerably larger number of users to be supported than the number 
permitted by the brute force IC scheme. Comparing Figures 10a, b with 
Figures 9a, b we observe that the optimum value of Pi obtained at a given 
value of degradation D is different for the two cases. 


Here again one can improve performance by going to a multiple stage 
scheme. For two stages of cancellation with linear tentative decisions as 
illustrated in Fig. 1 combined with Fig. 2b, it can be shown that the 
normalized output of the I&D corresponding to the final decisions of user 1 is 
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(by analogy with (37)) 
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( 66 ) 


Letting 


A, = P2 + A (1 - P2 ). Aj = AP2 (^ 7 ) 

A3 = a(i-a)+a(i-A2). a, = (i-aXi-A2) 

then the mean of the desired signal term (conditioned on Ao=D/ 
variances of the self noise, the user interference, and the equivalent thermal 
noise are respectively given by 
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Note that (69)-(71) reduce, respectively, to (38M40) when Pi = P2 = 1 since in 
this case, A, = Aj = 1 and P^ = P^=0. The corresponding degradation factor can 
be expressed in the form 
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The two stage linear partial interference cancellation scheme has also been 
computer simulated and the values of and p^ jointly optimized to provide 
for a minimum degradation. The performance results are illustrated in Fig. 
11 along with the corresponding results from Fig. 10a for a single stage 
scheme. 


Finally, two and three stage nonlinear (hard tentative decisions) partial 
interference cancellation schemes were computer simulated with optimized 
values of the partial cancellation paraimeters. The performance results are 
illustrated in Fig. 12 along with the corresponding simulation results for a 
single stage scheme. We observe that, for the parameters considered and 
uncoded BPSK users, a three stage nonlinear partial IC scheme allows as 
many as 80 users with a degradation of only 1 dB as compared to 9 users in a 
conventional CDMA system with the same degradation. This ideally 
represents an almost nine-fold increase in the user capacity of the system. 
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5.0 Parallel Interference Cancellation for CDMA with Coded Users 

Until now we have considered interference cancellation for CDMA systems 
in which the users transmit uncoded BPSK modulation. The IC techniques 
discussed in the previous sections are equally applicable to CDMA systems 
whose users employ error correction coding. To illustrate this applicability, 
we shall consider the performance of a CDMA system in which the users 
transmit identical convolutionally coded BPSK modulations. 

When used with a convolutionally coded modulation, there essentially 
exist two options for handling the "tentative" decisions made by each user in 
the presence of the interference. In one instance, these decisions can be made 
ignoring the fact that the each user’s data is encoded (that is, treating it as if it 
was an uncoded symbol stream). These symbol-by-symbol decisions would 
then be respread and remodulated directly on the carrier of each user. This is 
the simplest of the two options in that the only delay involved is the time 
(one coded symbol interval) to make these decisions. The second option is to 
treat the incoming data stream per user as coded s)anbol streams (as they truly 
are) and make the tentative decisions using a convolutional decoder (soft or 
hard decision) with the resulting bit stream being reencoded before 
respreading and remodulation on the carrier. This is the better (from a 
performance standpoint) of the two options but now involves a delay on the 
order of the decoder buffer size as well as a decoder/reencoder pair for each 
user at each IC stage. In view of this additional complexity, we shall consider 
only the first of the two options. 

The performance analysis of conventional (no IC) CDMA with 
convolutionally coded users was given in [18]. In this work, a constraint 
length K = 7, rate r = l/2 convolutional code with Viterbi decoding and 
eight-level soft decisions was assumed as the error correction 
coding /decoding scheme employed by each user. A summary of the 
performance of this code applied to PSK modulation in the presence of 
AWGN only , i.e., no multiuser interference, was given there in Table I. In 
particular, at a bit error probability Pt(E) = 10"*, the required coded bit energy- 
to-noise ratio is EJN^ = A.5 dB (as opposed to EJN^ = 9.6 dB for uncoded 
data) which results in a coded symbol energy-to-noise ratio 
Since the degradation factor only depends on the additional SNR required to 
overcome the multiuser interference and not on the specifics of the 
coding /decoding process, we can use the results of Section 3.1.1 as if we were 
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dealing with an uncoded system whose required bit energy-to-noise ratio is 
equal to the required symbol energy-to-noise ratio in the coded system. This 
is precisely what was done in [18] to assess the performance of conventional 
CDMA for coded users. In particular, for the above convolutionally encoded 
BPSK system, we have analogous to (16) 

(74) 

{N'oJi No + (M-1)5T; l + (M-l)77'V'‘(£./No)i, 
where r\ = T^IT^ is still the ratio of the chip rate to the information bit rate as 
previously defined and r = TJT, is the code rate with T, denoting the coded 
symbol duration. Since, however, the bit energy and coded symbol energy are 
related by £, = r£^, then (74) can be rewritten as 

(r£^/jVo)^ (75a) 

Jj 1 + (M - l)T?”V-'(r£^ /No), 

or 


( it.] (75b) 

UoJ, l + 

which is identical to (16). Hence, the degradation factor given m (17) is still 
applicable to the coded case. Since, however, in the coded case it is the coded 
symbols that are spread by the PN modulation, we shall in the numerical 
evaluations choose a value of tj equal to Umes the value used for the 
uncoded numerical results, e.g., for the rate 1/2 code, we shall use r] = 200, m 
order to maintain the analogy with the imcoded system. Once again, this is 
what was done in [18]. 

Insofar as the analysis for CDMA with IC is concerned, we can make the 
same analogy between the uncoded and coded cases as was done above since 
here again the IC process results from symbol-by-symbol decisions at each 
stage. Hence, as long as the Gaussian assumptions previously made for 
modelling the interference are valid, we can directly apply the imcoded user 
analytical results derived in the previous sections to tiie coded case. As far as 
simulation results are concerned, this allows using an encoder and decoder 
for only the desired user while treating the remainder of the users (the 
interferers) as imcoded. 

Fig. 13 is a plot (obtained by computer simulation) of the degradation 
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factor D versus the number of users M for a single stage IC scheme with 
hard tentative decisions. The value of p, has been optimized to yield the best 
performance (minimum D) and, in accordance with the above discussion, 
the values of /*(£) and tj have been chosen analogous to those in Fig. 4. 

6.0 Conclusions 

The inclusion of multistage parallel interference cancellation in a CDMA 
receiver can significantly improve its performance relative to that of a 
conventional CDMA receiver where no interference cancellation is 
attempted. A partial interference cancellation philosophy, in which the 
amount of interference cancelled is related to the fidelity of the tentative 
decisions involved in forming the interference estimate, is in general 
superior to a brute force philosophy of entirely cancelling the interference at 
each stage. Using a null zone type of device for making the tentative 
decisions at the various stages of the cancellation process is superior to using 
either a hard limiter or linear device. The linear device on the other hand 
has the advantage that the receiver implementation does not require 
knowledge of the user powers nor does it need carrier synchronization at the 
various stages. The latter implies that the final data decisions can be 
performed with a differential (rather than a coherent) detector. The 
technique is equally applicable to coded as well as uncoded modulations. 
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Figure 3 Tentative Decision Devices 
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(d) Null Zone Device 


Figure 3 Tentative Decision Devices (continued) 
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Fig. 8a. A Comparison of the Degradation Factors for One, Two, and Three Stage Brute Force Linear Interference 
Cancellation - Equal Power Users, /*j(£) = 10 
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Fig. 12. A Comparison of the Degradation Factors for One, Two, and Three Stage Nonlinear Interference 
Cancellation - Equal Power Users 
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Appendix A 


Analytical Models for Simulation of Multiuser 
Communication Receivers 

Consider the received signal of Eq. (1) corresponding to the transmission of 
M BPSK users over an AWGN channel in complex baseband form. For 
convenience this signal is rewritten here as 

K» =' + "(O = i^m, «)?«,«)«* +<■«) (*■■') 

i-1 i*l 

where ni(0 = f,a^i -kT,) is the data modulation of the i* user with unit 
power rectangular pulse shape p(t) of duration 7; (the bit time interval) and 

binary (±1) data bit in the transmission interval *:7; ^ t < (* + l)r„ n(t) 

is a complex white Gaussian noise process with power spectral density No 

watts/Hz, and PN,(t) is the user's ±1 PN code waveform. In the zero^ 
transmission interval 0 ^ t ^ 7],, Eq. (A.1) becomes 


r{t) = '^4s,a,oPN.{t)e’*‘+n{t) 

tsl 


(A.2) 


In deriving the decision metrics associated with the optimum and 
interference canceUation (IC) receivers, we shaU be interested in the 
projection of (A.2) on the M PN code waveforms treated as a set of basis 
fimctions. As such, we define the normalized received vector^ r with 
components 

r* r(r)PN,.(f)dr; i = 1,2,. ..,M (^'3) 

and the normalized AWGN vector n with components 

n. = ^\\{t)PN,it)dv, i = l,2 Af (A-4) 

The noise components (nj) are complex zero mean Gaussian random 


^All vectors defined will be assumed to be column vectors. 
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variables with variance = = and conditional 

crosscorrelation 


PNiit)PNjit)dt=N^Y^l ij = 1.2... .,M 

h 

where 




M 


(A.5) 


(A.6) 


are the components of the normalized PN crosscorrelation (user interference 
noise) matrix F . 


To proceed further, we model the i^ PN waveform over the zero^ bit 
interval 0 ^ r ^ T* as having been generated by a random binary (+1) data 
sequence {c^} of rate l/T^ (i.e., rj = T^/T^ PN code chips per data bit). Thus, we 
can write 


PiV.(0 = |^CaP(r-*7;) (A.7) 

where p(t)is agedn a unit power rectangular pulse shape now of duration T^. 
Defining the user MxN^ PN code matrix C = |c^ J whose row is the 
users code sequence in the interval Q^t^T^),i.e., C[ = (c„.c„.....c*^ ], then 
using (A.6). the PN crosscorrelation of (A.6) can be expressed as 


Yi/ ~ ^jk~ f~f *“ 1*2.....A/ (A.8) 

V *=i *-i V” 

and hence the user interference matrix F becomes 

r = icC^ = CC’' (A.9) 

where C = J is the normalized user PN code matrix. Note that the 
diagonal elements of F all have value imity. 

The noise components of (A.4) can also be expressed in terms of the user 
PN code matrix. Using (A.7) we have 





The noise components {N,] are independent zero mean complex Gaussian 
random variables with variance Alternately, 

defining the normalized noise vector C with zero mean, variance two {unit 

variance per dimension) Gaussian components = *en from 

(A.IO) the normalized noise vector n can be expressed as 


n = ^C 

Finally, defining the zero* bit interval data vector a. with components 
a ; / = 1,2,... ,M , the complex carrier phase diagonal matrix 


(A.11) 


o= 


0 


0 


gM 0 

0 • 


0 0 


0 

0 




and the scalar SNR diagonal matrix 


(A.12) 


B = 


Nn 


0 


0 


0 


Mi 

I N, 

0 


0 


’i>2 


0 


0 

0 


2^ 

iVn 


bM 


(A.13) 


u 

with E.MiT'b, we can write (A.l) in the normalized vector form 

(A.14) 

y = CC"^^Bao + CC 

where y=rl4NJ2 is the normalized received vector. For the case where all 
users have the same power, then (A.14) simplifies to 
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Appendix B 

Decision Metrics for Multiuser Communication 

Receivers 


/. Optimum (Maximum-Likelihood) Decision Metric — 1 Bit Observation Interval 

In this section, we derive the optimum decision metric for the multiuser 
signal of Eq. (1) based on an observation of the received signal r(r) over a single 

bit interval assumed, without loss in generality, to be the zero^^ transmission 
interval 0 ^ ^ In particular, we derive the metric for making a decision on 

the data vector ag = { 0 , 0 , 0 ^,... ,auof corresponding to the data bits transmitted 
by the M users in the interval 0^t<T^. This metric is a special case of that 
derived by Verdu [1] who considered maximum-likelihood sequence 
estimation of multiuser signals corresponding to an observation interval 
greater than a single bit in duration. Although, one can immediately reduce 
the results in [1] to the special case of interest here, our purpose is to obtain the 
results in the form of the models and notation set up in Appendix A and also 
to use them as a motivation for deriving the interference cancellation (IC) 
metrics that follow in subsequent sections. With this in mind, we shall 
independently derive the needed results here in this section. 

Using the well-known maximum-likelihood (ML) approach, the optimum 
decision rule for ag based on an observation of r(t) over the interval 0 ^ t ^ 7], is 
to choose the particvilar sequence minimizes the squared 

Euclidean distance 




i«l 


isl 


dt 


= [yiOfdt - 2Re| r(t)^4s,a,^PN.Xty-'*‘dt 

I i=l 


(B.l) 


Ignoring the first term in (B.l) since it does not depend on 

dividing the remaining two terms by two, we get the equivalent decision, 

namely, choose the particular sequence maximizes 
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(B.2) 


Q = Re|J^‘ r{t)^^4s,a,,PN,{t)e-^dt^ 

-|r[t ]<* 

Rewriting (B.2) in terms of the notation introduced in Appendix A, we get 


I i*l ^ J ^ 1=1 1=1 

= ^Wi |^W' 4 --iS 

2[ lr"V ^0 J 2trtfViV„ A^o 


(B.3) 


where as before £4, =5,T4 is the user's energy per bit. Alternately, in vector 
notation, we wish to maximize^ 


Q = Rc{a jBO'y} - ^ajBOD^'Bao 
= RejajB^Oy - |«>r«&*Bao]| 


(B.4) 


where the matrices ® and B are defined in (A.12) and (A.13), respectively. 
Because of the symmetry of the F matrix [see (A.6) and (A.9)], we can rewrite 
(B.4) in the simpler form 


C2 ^e{ajB«I>*y} - -ajB<I»*riI»Bae 

2 

= Re|ajBO*[y-ir«I»Bao]| 


(B.5) 


Again for the case where all users have the same power, maximizing (B.5) is 
equivalent to maximizing 


f2 = Rc 


aJO* 


y-- ^C'C’’«I»ao 
^ 2\ ” 


(B.6) 


II. Interference Cancellation (IC) Decision Metric Based Only on Observation of the 
Received Signal in a Single Bit Interval 

Suppose now that instead of making a joint decision on the M data bits 


^We shall continue to use the notation Q for all normalized forms of the decision metric. 
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transmitted in the zero*^ transmission interval, we instead choose to 
individually decide on each user's data bit in this same interval. The 
motivation for doing this is to reduce the complexity of the metric derived in 
the previous section. Clearly, in deriving such an ML metric for any one user, 
one would theoretically have to exact knowledge of the M - 1 data bits 
' corresjK>nding to all the other users. Since indeed this information is 
unknown, the above theoretical assumption is practically invalid. 

Nevertheless, let us proceed to find the appropriate ML metric under this 
assumption after which we shall modify it to allow for a practical realization in 
terms of >an IC decision algorithm. For simplicity of notation, we shall assume 
that user #1 is the user of interest, i.e., we seek an ML metric for deciding on 
data bit o,o- 

If once again we were to start by trying to minimize (B.l), assuming now 
that the sequence (o 2 o>%)>- ^ known a priori, then one would again 

arrive at the desire to maximize (B.3) with the terms that do not involve Ojo 
eliminated. Eliminating such terms from (B.3) gives the simplified decision 
rule: Choose the particular Ojo that maximizes 


1 . 


fiK 


u 


i2 = Re 
Note that 

theoreticcJly requires exact knowledge of the data sequence 
Before proceeding we rewrite (B.8) in the vector form 

X2= -(r-I),OBao]} = Rc[a,^-^[y, -(C’C^ -l) OBao]} (B.9) 

where I is an MxM identity matrix emd the "1" subscript on F— I and C C ^ - 1 
denotes the first row of these matrices. Since for a single user communication 
system, the optimum ML decision metric for 0,0 would be to choose the 
particular that maximizes 




-/♦l 


u 


i=2 


2£fc ;f 

-rr^iorue 

Nr, 


(B.8) 


as predicted above, a decision on based on maximizing (B.8) 


f2 = Re{oioe’'^^>'i} 


(B.IO) 



/ 


B-4 

then the subtractive term C 4»Ba in m 

r^rjr X . 

«*„ would analogou^^^^^ ‘^e decision on 

X2-Re{a e ^^[y^_(r_i)^OBao]}==Re[a^e-M[^^_(Cf^^ jj 

dateWt§Xo; / knowledge of the M - 1 

Since a^g in (B.9) is binary (+1) valued th#»n fho ^ • • 

can be formulated as a compaion oTdl' ^ 

equivalently, ^ threshold, or 

A.=agn{Re{e-M[,,-(Ccr-,),«^m 

where a,, denotes the decision on a, and "sen” h . v 
An analogous decision statisHc bas^ on (B llTw ^ 

user data bits in the Wh J„ ’ other 

deci^ons into a vector I wr^Tr^ ^^td^rr ^ ^ 

5. =5p'{Ra{'I''[y-(CC’'-I)®BaJU 

:r;xxx“crzr:ci'TX“: " t "■ 

reaUty,onlyaifectedbyftevectora.=(a,“.^“'r^ “”P°"7 “ 

ahtt^d^t ~7me 

deciding on a„. OneZ*!>d of » *®°«"«lge of a. m so far as 

acratm, (in time) approach as folloT’^^ to “ 

ueratton we replace a, on the right hand IHf (B 3) h 
namely, a,(0), obtained from tht soluHon , ,n ^ “f it, 

problem. In particular, 
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ao(0) = sgn{«]>*y} 

where we have introduced the parenthetical notation "(fc-1)" to denote the 

values of the corresponding variables at the beginning of the stage of 
iteration. (In terms of the present discussion, we have k = 1.) We refer to 
ao(0) as a "tentative decision" vector. 2 Note that each component of ao(0) 
requires no knowledge of the remaining user data bits and as such represents a 
decision made in the presence of the total interference from all other users. 
Needless to say these decisions will be poor relative to our final goal and it is 
for that reason that we refer to them as tentative decisions. Using ao(0) as a 
replacement for a„ on the right hand side of (B.13) we obtain after one stage of 

iteration the decision statistic^ 

ao(l) = sgn|Re{«&*[y - (C'C - l)«I>Bao(0)]}| (B15) 

Continuing this iteration for more than one stage will further improve our 
final decision on ao since now the tentative decision vector for the second stage 
will be ao(l) which has the benefit of one stage of interference canceUation. 

Thus, after k stages of iteration we obtain the final decision^ 

5„(t) = sgn{Re{*-[y - (C'C’- - l)4.Ba„(t - 1)]}} (B16) 

It is important to understand that the iterative decision statistic described by 
(B.16) or any variation thereof is not optimum in the ML sense and thus its 
performance will be inferior relative to that corresponding to the ML approach 
taken in Section I of this appendix. Nonetheless, as will be demonstrated in the 
main text, by proper design, iterative IC receivers based on the notions leading 
up to (B.16) can be made to yield performance comparable with the optimum 
ML receiver as suggested by Verdu but with considerably reduced implementation 
complexity. 

2Shortly, we shall modify this "tentative decision” to allow for a nmre general form of 
decision device. 

3ln reality, y here corresponds to the normalized received vector in tire zero* transmission 
interval delayed by one bit interval. Since, however, a delay of one or more bit intervals does 
not affect the evaluation of y (as given by (A.15)), we shall continue to denote by y this 
normalized received vector at the fcth stage of iteration, i.e., y delayed by k bit intervals. 
^Again, we shall soon allow for modifying the tentative decisions made at all stages of 
interference cancellation prior to the final decision in terms of a more general form of decision 
device. 
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Returning to the notion of the tentative decisions made at each stage of 
interference cancellation/ one should not conclude that a hard decision device 
wiU necessarily yield the best ultimate (after say K stages of IC cancellation) 
decision on Eq. Indeed, it will be shown in the main text that other forms of 
zero memory nonlinearity® for a tentative decision device can produce 
improved performance relative to the signum fimction. In fact, even a linear 
tentative decision has certain merits as will be discussed in the main text. In 
view of this, we shall generalize our iterative IC cancellation decision statistic 
to allow for an arbitrary tentative decision function which we shall denote as 
/(jc). Thus, in this more general context, (B.16) becomes^ 

ao(ik) = /(Re{<l^‘[y - (C'C"' - - l)]})=/(ao(*)) (B-17) 

or for the case where all users have equal power 


ao(*) = / 


1 

r r 

Re\ 

r 

1 

[ L 


y-. 


lEk 

N. 


(C'C'^-l)4»ao(*-l) 


=/(aoW) 


(B.18) 


11 uu uaci O ^ ^ \ j 

iteration) the signum function should be used as the decision device. Also, 

ao(Jt) represents the vector of tentative soft decisions at the stage of 
iteration. The decision metric in (B.18) suggests the implementation of an IC 
receiver as illustrated in Figure 1 together with Figure 2a. We refer to this 
structure as a brute force interference canceller since at each stage of the iteration 
an attempt is made to completely cancel the multiuser interference. 


Ill Interference Cancellation (IC) Decision Metric Based on Joint Observation of the 
Received Signal and the Tentative Soft Decision in a Single Bit Interval 

The decision metric derived in Section II as described by (B.17) implies that 
at each stage (iteration) the receiver attempts to entirely cancel the full amount 
of multiuser interference. Since in the early stages of interference cancellation, 
the tentative decisions are less reliable than they are in later stages, it is not 
intuitively clear that the above philosophy of entirely cancelling the full 

®In tfie next section of this appendix where the IC metric is generalized still furtirer to 
include additional observatiorw, we shall derive a specific form of this nonlinearity based 
upon nonlinear estimation considerations for Egf^). 

^Here again we use the notation /(x) to denote d\e application of the function to each 
component of the vector x . 
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amount of interference at each iteration stage necessarily leads to the best 
suboptimum decision metric. Rather, a better philosophy is one which m the 
early stages cancels only a fraction of the multiuser interference with the 
amoimt being cancelled increasing as one continues to iterate toward the 
ultimate final data decisions, i.e., as the fideUty of the tentative deasions 
improves. To see how a metric motivated by such a philosophy can come 
about using ML-type considerations, we proceed as follows. 

Analogous to Section H, we shall take an iterative approach to a 
suboptimum decision statistic derived from now joint ML-type considerations. 
In particular, consider first the component of the normalized received signal 
vector y corresponding to user #1 which when multiplied by is obtamed 
from Eq. (A.15) as 


Taking the real part of (B.19) gives an expression of the form 


(B.19) 


l2E, 




(B.20) 


N, 


Here 7 denotes the real part of the multiuser interference experienced by user 
#1 due to the remaining M - 1 users and W. is a zero mean Gaussian random 
variable with variance e[W^} = <r?. Next, consider the component of the vector 
of tentative soft decisions at the fc^h stage corresponding to user #1. From 
(B.17) together with (B.19) we have 


2,,(* - 1) = Reje-'V, - 

= Re| 1 

tVo i=2 V ^0 


(B.21) 


= jlE -H rJ X &(a«, - a„(k - 


+ W, 


For the purpose of what follows we can model the middle term of (B.21), i.e., 
the residual (uncancelled) interference, as a zero mean Gaussian random 
variable which when combined with W, results in a zero mean Gaussian 
random variable W^fJt-l) with variance E[W^(k-\)] = <Tlf Note that the 
variance al, depends on the iteration stage k. In view of the previous 
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discussion, we rewrite (B.21) in the form 


= + (B.22) 

Clearly, and are correlated, i.e., E{W^W 2 (A:-l)}=PtCr,cr 2 t. As we shall 

see shortly, it is not necessary to be able to specifically evaluate a„ ^ 2 * Pt- 
Rather a specific combination of these parameters will be used to define a 
parameter p* which shall have significance in terms of the amount of 
interference for which cancellation is attempted at each stage of the receiver. 


Using (B.20) and (B.22), a ML estimate of Ojo based on joint observation of Y^ 
and is obtained by choosing the value of that maximizes the 

conditional probability density function (pdf) 




Ina, 




Pi 


rexp 


’Ik 




2E, 


N, 


^0 


+ CT 


5 ,o(*-l)- 


2E, 


Nr 


“^0 


2a,V^t(l-Pt) 





2£k ^ 

{- „ 

2E. ^ 


^^0 J 

|flio(*-l)-J 

— -^0 
N 

^^0 y 


pI) 


(B.23) 

Since from (B.19) the interference term /, is a function of the data sequence 
- which is unknown, we proceed as before by replacing 7j with 


M 




i>2 


(B.24) 


which is synonymous with replacing the unknown sequence (o 2 o>%)’ - ’^vo) 
with the stage estimate of it Substituting 

/,(*) for /, in (B.23) and simplifying yields 
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p(r„a,o(*-l)|/,(A:)) = Cexp 


= Cexp 



|2£.. 

N, _ 




. 



A (^)) + ~ (Yi - A (^) + ^0 (^ - 1)) 

(al - Pt<rig2*)(y| -AW) + K -P*<^i^2t)^o(^-l) 




(B.25) 


where the constant C includes terms that do not depend on a,Q. We now 
introduce a normalization so as to make the coefficients of —I^(k) and 
o,o(it-l) in (B.25) sum to unity. Letting^ 


A olk-Pk^xOik 1 ^ Oyk-PkOl^lk 

* erf + a\t - “ ^Pk^^i^ut 

then (B.25) can be written in the desired form 


(B.26) 


p(y,.a,o(* - l)|aio. A(^)) = <^exp^ 



af+al-2p,a,a^lt 

No 



•(l-p*)a,o(*-l)] 


(B.27) 


Since the natural logarithm is a monotonic function of its argument, then 
taking the natural logarithm of (B.27), we see that the ML value of Ojq at the 

fcth stage of iteration is simply given by 


a,o(*) = sgn|[p*(i; - /i(fc)) + (1 - Pt)5,o(^ - 1)]} 


= sgn 


Re 


Pk 


Sv ^0 


>J 


+ (l-P*)flio(*-l) 


(B.28) 


fVi 

In terms of the entire vector of decisions for all users in the zero*-^ 
transmission interval, we have analogous to (B.16) 


^Note that based on its definition, the parameter is not necessarily restricted to lie in 


the range 0 ^ p* ^ 1. However, if < min — then it can be shown that this 

V^ 2 * cr, J 

restriction is valid. This implies either (T, ^(Tj* or Cj* <<Tj. Intuition, however, would 
suggest the former. 
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Ho(ife) = sgn{Re{p**I»‘[y - (C'C^ - l)OBa„(/: - 1)]} + (l - p*)a„(^ - 1)} (B.29) 

In applying (B.29), the parameter p* is to be optimized at each stage of the 
iteration. The decision metric of (B.29) suggests® the implementation of an IC 
receiver as illustrated in Figure 1 together with Figure 2b. 

The iterative joint ML approach that leads to (B.29) produces a decision 
metric in which, the tentative decisions at each stage of the iteration are hard 
decisions. We now consider a modification of this approach which allows for 
other than hard tentative decisions. (Of course, the final decision on a,o will, 
for uncoded data, always be a hard decision.) In particular, we consider using 
for fl,o(ik) the nonlinear estimate £|a,o|l^,a,o(^-l),/i}which is given by 


- 1)./,} = il)p{%a,oik - 1 )./,)+ (-l)p(-l|r„a,„(i -1),/,) (B.30) 

with p(a,o|y,,5,o(/k-l),/,) the a posteriori probability of the user #l's bit in the 

zero^^ transmission interval given the observations and interference. Using 
Bayes' rule, this probability can be determined in terms of the conditional 
probability p(yj,5io(I:-l)|a,o,/,) of (B.23) as 


p(aio|y„aio(^-l).A) = 


yi(l)<7(l) + yl(-l)<7(-l) 


where 


(B.31) 


>'(flio)=p(l'i.5,o(*-l)|aio.A). ‘7(oio) = p(^o.A) <B-32) 

Evaluating (B.31) at flio = 1 and a^Q = -1 and substituting the results into (B.30) 
gives 


As before, evaluation of (B.33) ideally requires knowledge of tiie interference /j 
which in turn requires that user #1 have knowledge of the sequence of data bits 
corresponding to the other users, i.e., In the absence of this 


exact knowledge we do as before and replace /j by the estimate /,(A:) as given by 
(B.24). As such we redefine the fvmctions in (B.32) by 


®As in Section II, we replace the hard tentative decisions with an arbitrary zero memory 
nonlinearity /(x). 
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^(flio)=p(i'i-^o(*-l)|^o.A(*)). <7 (Oio) = /’(^o.AW) 

Because of the symmetry of the problem, i.e., the equiprobable properties of the 
data streams, we have that ^(1) = <?(— 1). Hence, (B.33) becomes 

- D.w] = 

Referring to (B.27) for the evaluation of have that 

y\(fl,o) = Cexp(fl,oaA) 

where 


^ A IE^\ £tt^LzlMi£iL 

* “VNoL cr^<^u(^-pl) J (B.37) 

- /, {k)) + (1 - p,)a,,{k - 1) 

Finally, substituting (B.36) into (B.35) we obtained the desired result as 


;i a\ C exp(-«,Aj _ , . 

° Cexp(atA4) + Cexp(-atAt) * * 


Z \ al + cl, 2p,<^a A ^ 
' 0 ^Ik Pt j J 


(B.38) 


Comparing (B.38) with (B.28) we see that the hard tentative decisions have been 
replaced by soft tentative decisions in the form of hyp>erbolic tangent functions. 
Furthermore, the slope of these fxmctions (which is proportional to a*) is 
another parameter to be optimized at each stage of the iteration. 
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Appendix C 

The Special Case of Linear Interference Cancellation 


We consider here the special case of the interference cancellation methods 
discussed in Sections II and lU of Appendix B when the tentative decision 
devices are all linear, i.e., the tentative decisions are infinite bit soft decisions. In 
particular, we shall be interested in comparing the limiting (as the number of 
iterations approaches infinity) performance of such schemes with other linear IC 
cancellation methods such as those considered by Verdu [4]. 

For the brute force IC scheme as described in Section n of Appendix B, if we 
linearize the tentative decision devices, then (B.18) becomes^ 





^ 1 — -v‘ 


y-(r-I)ao(A:-l) 


(C.l) 


= y-(r-l)a,(k-l) 

where y denotes further normalization of the normalized received vector y. 
This equation when iterated can be expressed in the form 


1-1 


In the limit as /—>«», (C.2) becomes 


y+(-l)'(r-I)'ao(^) 


(C.2) 


a„=[I + r-I]-'y =r-'y (C.3) 

assuming of course that the inverse correlation matrix exists (more about 
this constraint shortly.) It is straightforward to show that the estimator in (C.3) is 
also the minimum mean-square error (MSE) linear estimator of the zero^^ 
transmission interval data vector So, i.e., if the normalized observation y is 

given by [see (A.15)] y = -j^C'C'^ao + CC or equivalently y = ra„ + 



^For simplicity of the presentation, we shall assume the case where all users have equal 
power and also ignore the effects of the user carrier phases. As such, the Re( ) notation in (B.18) 
can be ignored since the normalized received vector y is now real. Also, under these assumptions, 

the tentative decision transfer function becomes / (x) = {J'2-E^ / ■ 
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en tor of the form a„ = Ay, the choice of A = min* • fn 21 
Thus, we observe that the brute force IC scheme with }. 

asymptotically (as the number of stages r decisions 

the mean-square sense) linear estimate of U.e “>" *>«> (in 

Since d.e final decision must he a hard decision, .hen (C.3) is repla^d h, 
ao=sgn{l^y} 

In [4], the authors recognize the fact that p-* m 

decision algorithm based on (C.4) where p-* T ^ ^ 

generalized inverses of the crosscorrelation matrivT^^-^ ^ member of the set of 

1 ( )( for finite/ even teheni-. may not e«s. we can 

nT^hlr^f ™ag4T» <f«Uions and a fhUte 

pproxunanon to Uie imear mse estimator of (C.3). 

once again hnearize fte^Ve fte^SsbTblmr'^ “ "" 

So W = 


(C.5) 




CC *)[p*ao(^-l) + (l-p^)ao(it-2)]j 

y - (r - I)[p^ao (/: - 1) + (1 _ p^)ao(/: - 2)] 


When iterated, (C 5) can ho o */ o( 2)J 

/ can be expressed in the form 


ao(^ + 2/) = 




(C.6) 


‘idependenl of the values of the pe's. Thus ri.^ ^ >»«>mes equal to (C.3) 
tentative decisions and a finite numbe t ^ scheme with linear 

approximation to the linear MSE estinLto/*?^ ^ another 

C.2) and (C.6) approximate Se^S' ‘“'^®® between how 

converge toward this asymptotic solution as / -!1 " “ for **’®’' 

/ i.e., tor any fixed finite /, 
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with proper choice of the pj^'s, the generalized IC scheme will offer the better 
p>erformance of the two. 
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Appendix D 

Mean and Variance of Normalized I&D Output 
for Single Stage Linear Interference Cancellation 


The real part of the normalized output Xj, of (32), namely, consists of 
three terms. The first is the desired signal component which has mean 
(conditioned on Oio = 1) 


{ /' M \ , I Af-1 


(D.l) 


and variance 
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(D.2) 

Since for purely random PN codes, we have £{y,'‘} = 3/77^ -2/tj (see 
Appendix E). then 




21 




2(M-i)r n 

I ^J. 


(D.3) 


The middle term in (32) represents the residual user interference after 
cancellation which has zero mean and variance (see Appendix E) 
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The equivalent thermal noise has zero mean and variance 
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Finally, then the mean of (conditioned on Ojo =1) is given by (D.l), i.e., 

= P.6) 

and the variance of is given by the sum of (D.3), (D.4) and (D.5), namely. 


.2 A^n 


<T„ = 


, M-1 , (M-1)(M-2)1 . E, ( M-1) 

1 ' 7 « _2 


ri 
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M + 2 + 


M^-5M + 2 


(D.7) 


Although in principle the above exact analysis can be applied to a multiple 
stage linear IC scheme, the difficulty associated with evaluating the 
expectation of higher order products of the yj/s makes this approach 
impractical. Thus, to allow for analytical evaluation in the multiple stage 
case, we shall need to apply some simplification with regard to the statistical 
characterization of the y,y's, in particular, for large values of rj it is reasonable 
to approximate them as independent zero mean Gaussian random variables 
each with variance 1/t]. The argiunent for justifying this assumption is the 
same as that given in the discussion following Eq. (12) of the main body of the 
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paper which relates to an IC scheme with nonlinear (hard-limiting) tentative 
decisions. Thus, while we do not need to invoke this Gaussian assumption 
for the single stage linear IC scheme, i.e., the results in (D.l) - (D.7) are exact 
for random PN codes assigned to the users, we shall nonetiieless recompile 
these results based on the above assumption to allow their similarity to be 
used as a justification for employing the latter in the multiple stage linear IC 

analysis. 

Since evaluation of (D.l) depends only on the variance r,u ‘Ws result 
remains unchanged. In going from (D.2) to (DS), the only difference lies in 
the evaluation of £{y,‘}. For y„ a Gaussian (rather than binomial) random 

variable, we would have £{r;} = 3 {£{r?,}f =3/n' and thus (D.3) would 
become 


M 
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= E, 


2(M-1) 
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In (D.4), the expectation of the fourth-order product of the y/s partitions into 
two expectations of second-order products of the y/s resulting m 
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k^i /#i 
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Similarly, for the independent Gaussian assumption on the 
becomes 


(D.4') 
y./s, (D.5) now 
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Finally, (D.6) remains unchanged and (D.7) now 
(D.4’) and (D.5') which gives 


reflects the sum of (D.3'), 
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Appendix E 

Various Moments of the User Crosscorrelations 

In this appendix we evaluate various moments of the user crosscorrelation 
y.. defined in (A.6) and expressed in terms of the user PN code sequences in 

(A.8), i.e., 

1 ^ (E.l) 

where we have made use of the fact that the number of chips per bit N, is 
equal to the spreading ratio t) = TJT^ defined in the main text. Then, 




riz: 
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£{nri.ri .} = - , 
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Appendix F 


Some Useful Closed Form Definite Integrals 


In this appendix, we tabulate some dosed form definite integrals that are 
particularly useful in evaluating the moments of (13). Many of these which deal 
with integrands containing a Gaussian probability density function (pdf) or 
Gaussian probability integral cannot be found in standard tables of integrals. 

Define the one- and two-dimensional zero mean Gaussian pdf's by 


Gi(x;<t,)= ^ 




(F.l) 


2na 




-exp<- 


<Tlx^ + aly^-2pa,<JyXy \ 

and the normalized Gaussian probabiUty integral [see Eq. (12)] 


G,(x,y;cr,,(T )= 
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Denote the hard limiter and null zone limiter functions, respectively, by 


fl, x>0 


nsgn X = s 


1, x>C 
0 , -C^x<C 
-1, x<-C 


= x<0’ 

Then, the following one- and two-dimensional integrals are appropriate; 
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One-dimensional integrals 

£exp{-A(l -t- Bxf}G,{x\a,)djc = ^ 1 + 2alAB ^ } 


(F.4) 



F-2 


f_xcxp{-Ax^ ± Bx}G,{x;a^)dx = expj— + ^ 

j^x^ exp{-Ax^ ± Bx}G,{x;a^)dx = 4-|-]exp|£|; ^ 


£xsgn(Ax + B)G^{x;(T^)dx = 
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Two-dimpnsinnal integrals 


£E exp |-[a(A + xf +b{B + }')^]|G2 (x;a,o)dxdy 

(oA^ +b^){l-p^) + (A^ + B^ -2pAB)abr\ 
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LLjcsgn(A + y)G2{x;a,,c^)dx = ^ ^ * exp^ 


2<Jt 
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